The hard to soft Pomeron transition in small x 
DIS data using optimal renormalization 

Clara Salas 

Instituto de Fisica Teorica UAM/CSIC and Universidad Autonoma de Madrid, 
E-28049 Madrid, Spain 

Abstract. We show that it is possible to describe the effective Pomeron intercept, determined 
from the HERA Deep Inelastic Scattering data at small values of Bjorken x, using next-to-leading 
order BFKL evolution together with coUinear improvements. To obtain a good description over the 
whole range of we use a non-Abelian physical renormalization scheme with BLM optimal scale, 
combined with a parametrization of the running coupUng in the infrared region. 



INTRODUCTION 

The description of Deep Inelastic Scattering (DIS) data for the structure function 
FiixiQ^) in different regions of Bjorken x and virtuality of the photon is one of 
the classical problems in perturbative QCD. The recent HERA combined results [1] for 
F2 cover a broad range of values of the photon virtuality and also reaches very small 
values of x, making this observable suitable to study the high energy limit, given when 
the center of mass energy of the system is much higher than any other scale involved or, 
in other words, the region of low Bjorken x. 

The aim of the present study is to analyze the structure function F2 and Pomeron in- 
tercept within perturbative QCD in this mentioned region of low values of x using high 
energy factorization [2], which convolutes the proton and photon impact factors with the 
gluon Green's function. For the latter we use the solution to the BFKL evolution equation 
at Next-to-Leading-Logarithmic (NLL) accuracy [3, 4] and introduce coUinear improve- 
ments, needed to deal with the leading coUinear singularities that are numerically large 
in this kinematic region [5, 6, 7, 8, 9]. We show how the use of a non-Abelian physical 
renormalization scheme with optimal scale setting^ allows for a good description of the 
Pomeron over the full range of Q^. 

A numerical analysis of the Pomeron intercept A and comparison with data is given, 
using the prediction of Regge theory that claims that in the high energy limit the 
asymptotic expression for F2 should grow with energy as (s/sq)^ = (l/x)^ . This 
analysis is presented in more detail in [12]. 



We use in this work MOM scheme with BLM optimal scale [10, 11] but there are other similar choices. 
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FIGURE 1. High energy factorization (left = a), proton impact factor (center = b) and model for the 
running coupling freezing in the IR (right = c). 



THEORETICAL SETUP 

High energy factorization allows to write the proton structure function as a convolution 
in transverse momentum space of a non perturbative object describing the proton (proton 
impact factor ^p) with the photon (photon impact factor 4>y), calculated using pertur- 
bation theory, together with a gluon Green's function /, linking both process-dependent 
components and accounting for the BFKL evolution, as shown in fig. la: 

F2{x,Q^) = (2^)4 y -kf y ^<i>,(k,)/(x,k«,k,)o^(k,) . 

We use for simplicity in our analysis the leading order photon impact factor as presented 
in [13]. Our choice for the proton impact factor is a Poissonian-like distribution given 
by (see Fig. lb) 

<I>p(k)=A,,(kVeo)'e-'^'/Qo 
Finally, the gluon Green's function is governed by the BFKL equation. Since in DIS 
^ 2o' ^^'^ '^^ written in a form consistent with the resummation of ajog (1/x) 
contributions: 

where K is the NNL BFKL kernel. The zeros of the denominator in the integrand 
generate in the limits 7— )► 0, 1 all-orders terms not compatible with DGLAP evolution [5, 
9]. The first of these pieces (^(a^)) is removed when the NLO correction to the 
BFKL kernel is taken into account but not the higher order ones, which remain and are 
numerically important. A scheme to eliminate these spurious contributions [5] consists 
of modifying the BFKL kernel by making a shift of the form K{y) — )■ K{y+ co/2). 

Once we include NLL corrections with the coUinear improvement as described in [9] 
we find the following expression for F2 in Mellin space: 

F2{x,Q^) - |~dv^-^(3+-)c2(v)cp(v)|l + a^og('l') J|;^o(^ + /v') (1) 



-log — - +M ;rcoth(;rv) -2;rtanh(;rv) -M2(v) 



where C2(v) and cp{v) are the photon and proton impact factors in v-space, respectively 
and M2 is a simple function of V defined in [12]. Note that in this equation we have 
exponentiated only the scale invariant terms of the kernel. The terms that break the 
invariance are due to the running of the coupling and appear as a differential operator in 
V acting on the impact factors [14]. 

We introduce the running of the coupling in a way that removes the jU dependence of 
eq. (2) by making the replacement 



This resummed coupling is consistent with the Landau one up to NLO accuracy. We see 
in our analysis that in order to have a good description of the Pomeron intercept over 
the full range 1 GeV^ < < 300 GeV^, we need to move to a renormalization scheme 
based on the existence of a possible IR fixed point. The pieces of the NLL BFKL kernel 
proportional to j8o are absorbed in a new definition of the running coupling so that all the 
vacuum polarization effects from the jSq function are resummed. We refer to the reader 
to [10, 11] for details. 

As a final ingredient, we use a simple parametrization of the running coupling intro- 
duced by Webber in [15] which is consistent with global data of power corrections to 
perturbative observables and with the usual running coupling with Landau pole up to 
NLO accuracy. Fig. Ic shows both models, being the solid curve the new one. 



NUMERICAL ANALYSIS & COMPARISON TO DIS DATA 

To obtain our theoretical results we have calculated the logarithmic derivative f^^^^ij^^ 
using Eq. (2) with some modifications. For the comparison with DIS data we chose 
the values Qq = 0.28GeV and 5 = 8.4, = 4 and A = 0.21 GeV (see fig. 2b). The 
experimental input has been derived from the combined analysis performed by HI and 
ZEUS in Ref. [1] with x < 10~^. In the results indicated with "Real cuts" we have 
calculated the effective intercept for F2 at a fixed Q^, averaging its values in a sample of 
X space consistent with the actual experimental cuts in x. To generate the continuous line 
with label "Smooth cut" we have used as boundaries in x space those shown in Fig. 2a, 
which correspond to an interpolation of the real experimental boundaries. Note that the 
difference between both approaches is very small. We would like to stress the accurate 
description of the combined HERA data in our approach, in particular at very low values 

ofe". 

It is possible to improve the quality of our fit by introducing subleading contributions 
such as threshold effects in the running of the coupling, heavy quark masses and higher 
order corrections to the photon impact factor which became recently available [16]. We 
are presently working on these improvements, together with a comparison to data not 
averaged over x, also including an analysis of Fj^. 




q'' (GeV^) q'' (GeV^) 



FIGURE 2. Experimental bounds in x (left = a) and evolution of A with (right = b). 
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